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In this paper quasi-stationary, two-and-a-half-dimensional magnetic reconnection is studied in the 
framework of incompressible resistive magnetohydrodynamics (MHD). A new theoretical approach 
for calculation of the reconnection rate is presented. This approach is based on local analytical 
derivations in a thin reconnection layer, and it is applicable to the case when resistivity is anoma- 
lous and is an arbitrary function of the electric current and the spatial coordinates. It is found that a 
quasi-stationary reconnection rate is fully determined by a particular functional form of the anoma- 
lous resistivity and by the local configuration of the magnetic field just outside the reconnection 
layer. It is also found that in the special case of constant resistivity reconnection is Sweet-Parker 
and not Petschek. 



PACS numbers: 52.30.Cv, 52.35.Vd, 52.65.Kj 



I. INTRODUCTION 



Magnetic reconnection is the physical process of break- 
ing and rearrangement of magnetic field lines, which 
changes the topology of the field. It is one of the most 
fundamental processes of plasma physics and is believed 
to be at the core of many dynamic phenomena in labo- 
ratory experiments and in cosmic space. Unfortunately, 
in spite of being so important, magnetic reconnection 
is still relatively poorly understood from the theoreti- 
cal point of view. The reason is that plasmas usually 
have very high temperatures and low densities. In such 
plasmas, the Spitzer resistivity is extremely small and 
magnetic fields are almost perfectly frozen into cosmic 
plasmas. As a result, simple theoretical models, such as 
the Sweet-Parker reconnection model [1, 2] predict that 
the magnetic reconnection processes should be extremely 
slow and insignificant throughout the universe. On the 
other hand, astrophysical observations show magnetic re- 
connection tends to be fast and is likely to be the primary 
driver of many highly energetic cosmic processes, such as 
solar flares and geomagnetic storms. This contradiction 
between theoretical estimates and astrophysical observa- 
tions triggered multiple attempts to build a theoretical 
model of fast magnetic reconnection. 

First, in 1964 Petschek proposed a fast reconnection 
model [3] , in which fast reconnection is achieved by intro- 
ducing switch-off magnetohydrodynamic (MHD) shocks 
attached to the ends of the reconnection layer in the 
downstream regions and by choosing the reconnection 
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layer length to be equal to its minimal possible value 
under the condition of no significant disruption to the 
plasma flow. However, later numerical simulations and 
theoretical derivations did not confirm the Petschek theo- 
retical picture for the geometry of a reconnection layer [4- 
6]. Second, numerical simulations studies of anoma- 
lous magnetic reconnection, for which resistivity is en- 
hanced locally in the reconnection layer, were pioneered 
by Ugai and Tsuda [7, 8], by Hayashi and Sato [9, 10], 
and by Scholer [11]. Third, Lazarian and Vishniac pro- 
posed that fast reconnection can occur in turbulent plas- 
mas [12], although the back-reaction of magnetic fields 
can slow down reconnection in this case [13]. Finally, 
recently there have been number of attempts to explain 
the fast magnetic reconnection by considering non-MHD 
effects [14-21]. 

Most previous theoretical and numerical studies con- 
centrated on reconnection processes in two-dimensions 
or in "two-and-a-half-dimensions" . The later is the term 
used for a problem in which physical scalars and all three 
components of physical vectors depend only on two spa- 
tial coordinates (e.g. x and y) and are independent of the 
third coordinate (z). 

In this paper we consider two-and-a-half-dimcnsional 
magnetic reconnection with anomalous resistivity in 
the classical Sweet-Parker-Petschek reconnection layer, 
which is shown in the left plot in Fig. 1. The recon- 
nection layer is in the x-y plane with the y-axis being 
along the layer and the x-axis being perpendicular to the 
layer. The length of the layer is equal to 2L'. Note 
that L' is approximately equal to or smaller than the 
global magnetic field scale, which we denote as L. The 
thickness of the classical reconnection layer, 2<5 Q , is much 
smaller than its length, i.e. 2S <C 2L'. The classical 
Sweet-Parker-Petschek reconnection layer is assumed to 
possess a point symmetry with respect to its geometric 
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FIG. 1: (Color online) The geometrical configuration of the 
classical Sweet-Parker-Petschek reconnection layer is shown in 
the left plot. Petschek shocks exist only if the reconnection is 
considerably faster than the Sweet-Parker reconnection rate. 
The right plot shows an enlarged picture of the central layer 
region. 



center point O and reflection symmetries with respect to 
the axes x and y (refer to Fig. 1). Thus, for example, the 
x- and y-components of the plasma velocities V and of 
the magnetic field B have the following simple symme- 
tries: V x (±x,tv) = ±V x (x,y), V y (±x,Ty) = TV y (x,y), 
B x (±x,Ty) = TB x (x,y) and B y (±x,Ty) = ±B y (x,y). 
There might be a pair of Petschek shocks attached to 
each of the two reconnection layer ends in the down- 
stream regions (see Fig. 1). Because of the MHD jump 
conditions on the Petschek shocks, the presence of these 
shocks requires the presence of a significant perpendicular 
magnetic field B x at the reconnection layer ends [5, 22]. 
The plasma outflow velocity from the reconnection layer 
is approximately equal to the Alfven velocity Va (if the 
plasma viscosity is not very large). The plasma inflow 
velocity Vr outside of the reconnection layer, at point 
M in Fig. 1, is much smaller than the outflow velocity, 
Vr <C Va- Finally, the magnetic field outside the recon- 
nection layer is mostly in the direction of the layer (i.e. in 
the y-axis direction). 

The problem of quasi-stationary anomalous magnetic 
reconnection in the classical Sweet-Parker-Petschek re- 
connection layer was recently theoretically addressed by 
Kulsrud [5] for the special case of zero guide field (B z = 
0), zero plasma viscosity and anomalous resistivity that 
is a piecewise linear function of the electric current [see 
Eq. (32)]. There are two major results of the Kulsrud 
anomalous reconnection model. First, in the case of the 
Petschek geometry of the reconnection layer and a con- 
stant resistivity, one has to calculate the half-length of 
the reconnection layer V from the MHD equations and 
the jump conditions on the Petschek shocks, instead of 



treating L' as a free parameter (as Petschek did) . When 
the layer half-length L' is calculated correctly, it turns 
out to be approximately equal to the global magnetic field 
scale, L' » L. In this case the Petschek reconnection re- 
duces to the slow Sweet-Parker reconnection [5, 22]. This 
theoretical result is in agreement with the results of nu- 
merical simulations of two-dimensional reconnection with 
constant resistivity [4, 6]. The second major result of the 
Kulsrud reconnection model is that in the case when re- 
sistivity is anomalous and enhanced (e.g. by plasma in- 
stabilities), the reconnection rate becomes considerably 
faster than the Sweet-Parker rate. 

In this paper we develop and use a new theoretical ap- 
proach for calculation of the reconnection rate in the case 
of anomalous resistivity. This approach is based on ap- 
plication of the MHD equations in a small region, which 
is localized at the geometric center of a thin reconnec- 
tion layer (point O in Fig. 1) and has a size of order of 
the layer half-thickness 5 . It turns out that by using lo- 
cal analytical calculations in a thin reconnection layer, we 
can derive an accurate and rather precise estimate for the 
reconnection rate. In particular, we find the interesting 
and important result that a quasi-stationary reconnec- 
tion rate is fully determined by the anomalous resistivity 
function and by the magnetic field configuration just out- 
side the reconnection layer (at point M in Fig. 1). The 
underlying physical foundations of our analytical deriva- 
tions are, of course, the same as those previously used 
by others [1, 2, 5]. However, our theoretical and com- 
putational approach is somewhat different from the con- 
ventional approach to the reconnection problem, and we 
explain the difference in the next section. There are three 
main benefits of our novel theoretical approach to mag- 
netic reconnection. First, our approach allows us to ex- 
tend the results for anomalous reconnection obtained by 
Kulsrud in 2001 [5] to a general anomalous reconnection 
case, when the guide field and plasma viscosity are arbi- 
trary and anomalous resistivity r\ is an arbitrary function 
of the electric current and the two spatial coordinates. 
Second, it gives a new important insight into the recon- 
nection problem, such as dependence of the reconnection 
rate on magnetic field configuration just outside the re- 
connection layer. Third, our approach, based on local 
calculations, is applicable to cases when there is no well- 
defined global magnetic field structure, such as the case 
of multiple current sheets in a turbulent plasma. In our 
calculations we make only a few assumptions, which are 
described in detail in the beginning of Sec. III. 

This paper is organized as follows. Because our deriva- 
tions are rather complicated in the general case of re- 
connection with anomalous resistivity, we find it useful 
and instructive to consider the Sweet-Parker reconnec- 
tion first and to compare our theoretical approach to the 
classical Sweet-Parker calculations in the next section. 
In Sec. Ill, we derive our general equations for magnetic 
reconnection with anomalous resistivity, including the 
equation for the reconnection rate. In Sec. IV we consider 
and analyze different special cases of magnetic reconnec- 
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tion, in which our equations simplify and become easier 
for analysis and comparison to the previous theoretical 
and simulation results. In Sec. V we present the results 
of our numerical simulations of unforced anomalous mag- 
netic reconnection. These simulations are intended for a 
demonstration of the predictions of our theoretical recon- 
nection model for the special case of the piecewise linear 
resistivity function that was considered by Kulsrud [5]. 
Finally in Sec. VI we give our conclusions and discuss our 
results. Some derivations are given in the appendices of 
the paper. 



II. THE SWEET-PARKER MODEL OF 
MAGNETIC RECONNECTION 

For simplicity and brevity hereafter in the paper for all 
electromagnetic variables we use physical units in which 
the speed of light and four times n are replaced unity, 
c = 1 and 47r = 1. To rewrite our equations in the stan- 
dard CGS units, one needs to make the following sub- 
stitutions for electromagnetic variables in the equations: 
magnetic field B — > B/V^r, electric field E — > cE/V^r, 
electric current j — > (\/47r/c)j, resistivity 77 — s- 77 (does not 
change), magnetic field vector potential A — > A/V47T. 

In this section we consider Sweet-Parker reconnec- 
tion [1, 2]. We assume that resistivity is constant, 
77 = const = r] , the plasma viscosity is zero, the guide 
field is zero (B z = 0), and the geometry of the recon- 
nection layer is the classical Sweet-Parker geometry with 
the layer half-thickness 5 <C L and the layer half-length 
V = L, as shown in Fig. 1. The purpose of this sec- 
tion is to introduce and explain our new theoretical ap- 
proach to calculation of the reconnection rate and to com- 
pare it with the classical Sweet-Parker calculations. For 
this purpose, we first present the classical, conventional 
derivation of the Sweet-Parker formula for the reconnec- 
tion rate, and afterward we present our new derivation of 
the formula and explain the difference between the two 
derivations. 

In the case of a quasi-stationary reconnection the mag- 
netic field in the reconnection region changes slowly in 
time, <9B/9i ps 0. Therefore, in the two- and- a- half- 
dimensional geometry (d/dz = 0) the Faraday's law 
equation — V x E = dH/dt ps results in the z- 
componcnt of the electric field being constant in the re- 
connection region, VE Z = and E z = E z (t) is a function 
of time only. On the other hand, in two- and- a- half dimen- 
sions the x- and y-components of Faraday's law equation 
dB/dt = - V x E = -V x (rjj - [V x B]) reduce to the 
following equation for the z-component of the magnetic 
vector potential A: 

-E z (t) = 8A z /dt = -(V • V)A Z + V V 2 A Z 

= V X By-V y B x -r, ]z . (1) 

(The resistivity r\ = const = i] is constant in the Sweet- 
Parker reconnection case, but equation (1) is general 



and is valid even if resistivity 77 is anomalous and non- 
constant.) 

Now, the left- and right-hand-sides of equation (1) 
are constant in space. Therefore, the right-hand-side 
of equation (1) is constant across the reconnection layer 
(i.e. along the x-axis). We equate its values at points O 
and M, which are on the x-axis and are shown in Fig. 1. 
As a result, we immediately obtain 

riojo = V R B m , (2) 

where we use the following notations: Point O is the 
geometric center of the reconnection layer, where the z- 
component of the current is j a = j z (x = 0, y = 0) and the 
plasma velocity is zero. Point M is a point on the x-axis 
just outside the reconnection layer, where the resistivity 
term can be neglected in equation (1), see Fig. 1. We 
also use the notations B y = B m and V x = —Vr for the 
y-component of the field and x-component of the plasma 
velocity at point M, and take the reconnection velocity 
Vr positive. Note that V y — B x = at point M because 
of the symmetry of the problem with respect to the x- 
axis. Next, we estimate current j a at the central point O 
of the reconnection layer as 

Jo « B m /5 a , (3) 

where S a is the reconnection layer half-thickness, and 
we use Ampere's law, j a = (dB y /dx) — (dB x /dy) ps 
(dB y /dx) ps B m /S , in which we drop the (dB x /dy) 
term because the reconnection layer is thin [36] . 

Next, consider the x- and y-componcnts of the equa- 
tion of plasma motion. We will see below that the 
Sweet-Parker reconnection is slow, Vr -C Va = B m j ' ■ s fp. 
Therefore, in the equation of plasma motion along the x- 
axis (i.e. across the reconnection layer) the inertial term 
can be neglected, and this equation becomes the force 
balance equation, (d/dx)(P + B 2 /2) = 0, resulting in 
P = P m + -Bm/2) where P and P m are the values of 
the plasma pressure at points O and M. Here we use the 
fact that the magnetic field is zero at the central point O 
because of the symmetry of the problem. As far as the 
equation of plasma motion along the y-axis (i.e along the 
layer) is concerned, the y-components of the magnetic 
tension force and the pressure gradient force are approx- 
imately equal in the Sweet-Parker reconnection case. In- 
deed, on the y-axis the tension force can be estimated as 
(B • V)B y = B X x (dBy/dx) ps {5 /L)B m x {B m /S ) = 
Bf n /L ps (P - P m )/L ps -dP/dy. Here we estimate the 
x-component of the field as B x ps (5 / L)B m because B x 
is produced by the rotation of the B y field component in 
the reconnection layer [5] , and in the Sweet-Parker model 
it is assumed that the downstream pressure is equal to 
the upstream pressure P m (sec [1, 2]). The pressure and 
magnetic tension forces accelerate plasma along the re- 
connection layer (i.e. along the y-axis) up to the down- 
stream velocity V out , which can be estimated from the 
energy conservation equation 

(l/2) P V 2 ut ps B 2 J2 V out *V A = B m Up. (4) 
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This equation means that the work done by the pressure 
and magnetic tension forces along the entire reconnection 
layer is equal to the kinetic energy of the plasma in the 
downstream regions. 

Finally, in the Sweet-Parker reconnection model the 
plasma is assumed to be incompressible. Therefore, the 
mass conservation condition for the entire reconnection 
layer results in 

LV R ^S V ouU (5) 

where V out , the velocity of plasma outflow in the down- 
stream regions, is given by equation (4). Using equa- 
tions (2)- (5), we obtain the formula for the Sweet-Parker 
reconnection velocity [1, 2], 

V R m Va^o/VaL) 1 ' 2 . (6) 

Equations (2)-(6) are the classical Sweet-Parker equa- 
tions. Note that equations (2), (3) and (6) are local, in 
the sense that they are written for a small region of space, 
which is localized at the geometric center of a thin recon- 
nection layer (point O in Fig. 1) and has a size of order 
of the layer half-thickness S . All physical quantities that 
enter these three equations are defined in this small re- 
gion of space. At the same time equations (4) and (5) are 
global, in the sense that they result from consideration of 
the entire reconnection layer and they include the plasma 
outflow velocity U ou t , which is a physical quantity in the 
downstream regions at the ends of the reconnection layer. 
In our new theoretical approach to the calculation of the 
reconnection rate we intend to use only local equations. 
Our intent and our derivations will be justified by the new 
and important results that we obtain in the next section 
of this paper and discuss in Sec. VI. At present, let us 
explain our theoretical approach for the simple case of 
Sweet-Parker reconnection that we consider in this sec- 
tion. 

In the derivation of our theoretical model for magnetic 
reconnection we keep equations (2) and (3) unchanged 
because these equations are local. However, we rewrite 
equations (4) and (5) because they are global. To rewrite 
these two global equations in a local form, we consider a 
point O = (x = 0, y = y) that is located on the y-axis 
in an infinitesimal vicinity of the reconnection layer cen- 
tral point O (see Fig. 1) and has an infinitcsimally small 
value of its y-coordinate y. Since y — > +0, along the 00 
interval we can use the up-to-the-first-order Taylor ex- 
pansions B x (0,y) = (3y, V y (0,y) = vy and j z (0,y) = j Q 
for the values of the perpendicular magnetic field B x , 
plasma velocity V y and z-component of the current j z . 
These expansions are along the y-axis, and, of course, 
(3 = (dB x /dy) and v = (dV y /dy) are the first-order 
partial derivatives of B x and V y at point O (note that 
(dJ z /dy) — 0). Now equation (4) for the plasma accel- 
eration along the y-axis can easily be rewritten in a local 
form at point O as 

\p{vyf « £ B x (0,y)j z (0,y)dy = £ I3yj dy 



= f3 3o f/2 pv 2 w 0j o , (7) 

where the left-hand-side of this equation is the plasma 
kinetic energy at point O. On the right-hand-side of 
equation (7) we keep only the magnetic tension force 
for plasma acceleration because, as we found above, in 
the Sweet-Parker reconnection case the y-components of 
the magnetic tension and pressure gradient forces are ap- 
proximately equal to each other (note that in our general 
calculations in the next section we will take all forces into 
account). Next, equation (5) for the mass conservation 
of an incompressible plasma can easily be rewritten in 
the local form inside the area OMMO shown in Fig. 1 as 
yV R rs S o V y (0,y) = S a vy. Thus, we have 

v « V R /S , where v = {dV y /dy) = -{dV x /dx) . (8) 

This equation can also be viewed as the first order Tay- 
lor expansion of V x along the x-axis (i.e. across the re- 
connection layer), Vr = — V x (5 o ,0) w — (dV x /dx) 5 = 
(dV y I 'dy) S = vS , where we use the plasma incompress- 
ibility condition dV x /dx + dV y /dy = 0. 

Note that our first and second equations are local and 
are exactly the same as equations (2) and (3) in the 
Sweet-Parker reconnection model. Our third and fourth 
equations (7) and (8) are also local and are the analogues 
of the two global Sweet-Parker equations (4) and (5). Us- 
ing local equations instead of global ones is the first major 
difference between the Sweet-Parker and our theoretical 
models. However, in our theoretical model we have an 
additional unknown parameter [3 = (dB x /dy) Q1 which 
does not directly enter the classical Sweet-Parker recon- 
nection model. As a result, we need one additional equa- 
tion to be able to calculate the reconnection rate under 
the framework of our local model. This additional equa- 
tion comes from the condition that the right-hand-side of 
equation (1) is constant not only across the reconnection 
layer, but also along the layer that is along the y-axis. 
This condition is not explicitly used in the Sweet-Parker 
model, but it is used in our model, and this is the second 
major difference between the two models. To derive the 
additional equation, used in our reconnection model, we 
differentiate the left- and right-hand-sides of equation (1) 
along the y-axis (i.e. along the reconnection layer). The 
first order partial derivatives d/dy are identically zero 
because of the symmetry of the problem with respect to 
the x-axis. The second order partial derivatives d 2 / dy 2 
of the left- and right-hand-sides of equation (1) result in 

= -2{dV y /dy) (dB x /dy) -i lo {d 2 ]z /dy 2 ) 
« -2vp + 2 VoJo /L 2 , (9) 

where we use the fact that V x = B y = on the y- 
axis, we use our definitions (3 = (dB x /dy) and v = 
(dV y /dy)o, and we estimate the second derivative of the 
current j z as (d 2 j z /dy 2 ) w —2j /L 2 . Now all five equa- 
tions (2), (3), (7)- (9) are local. Combining them together, 
we easily obtain the Sweet-Parker formula for the recon- 
nection velocity, given by equation (6), which is naturally 
also local, see Ref. [37]. 



■5 



The reader of this paper could question why we develop 
and suggest a new theoretical approach to the problem 
of quasi-stationary magnetic reconnection if we obtain 
the same results in the Sweet-Parker reconnection case. 
The answer is that our approach, based on local calcu- 
lations, allows us to calculate the reconnection rate in 
the case of anomalous resistivity and also provides an 
additional important understanding of the reconnection 
problem. Our derivations for anomalous reconnection are 
given in the next section and we discuss our results in 
Sec. VI. Note that our local-equations approach to the 
reconnection problem and the more conventional global- 
equations approach are the same from the point of view 
of the underlying physics. Indeed, it is well known that 
by using the Gauss-Ostrogradski theorem, most physical 
equations can be written in two equivalent forms, in the 
form of local differential equations and in the form of 
global integral equations. 



III. MAGNETIC RECONNECTION WITH 
ANOMALOUS RESISTIVITY 

In this section we study reconnection with anomalous 
resistivity and derive a simple and accurate estimate 
of the reconnection rate in the classical Sweet-Parker- 
Petschek two-and-a-half dimensional reconnection layer 
shown in Fig. 1. We consider resistivity to be a given 
arbitrary function of the z-component of the electric cur- 
rent and the two-dimensional coordinates, r] = rj(j z ,x, y), 
which has finite derivatives in y up to the second order 
and in x and j z up to the first order [38]. 

Let us list the assumptions that we make. First, we 
assume that the characteristic Lundquist number of the 
problem is large, which (by our definition) is equivalent 
to the assumption that resistivity is negligible outside the 
reconnection layer and the non-resistive MHD equations 
apply there. Second, we assume that the plasma flow 
is incompressible, divV = 0. In the limit of very high 
Lundquist numbers and slow reconnection rates the in- 
compressibility condition is a very good approximation in 
a reconnection layer even for compressible plasmas [23]. 
Third, we assume that the reconnection process is quasi- 
stationary. This can only be the case if the reconnection 
rate is small, ij j /V A B m w V R /V A < 1 [V A = B m /^/p 
and see Eq. (2)], and there are no plasma instabilities in 
the reconnection layer. Note that in our model the re- 
connection rate can still be much faster than the Sweet- 
Parker rate. Fourth, we assume that the reconnection 
layer is thin, S /L' <C 1. If plasma kinematic viscosity is 
small (in comparison with resistivity), then the plasma 
outflow velocity in the downstream regions is equal to 
the Alfven velocity V A , we have plasma mass conserva- 
tion condition S /L' Vr/V a , and this assumption of a 
thin reconnection layer is fully equivalent to the previ- 
ous assumption of a small reconnection rate. However, 
if plasma is very viscous, then the plasma outflow ve- 
locity is smaller than V A and assumption 8 /L' <C 1 is 



stronger than assumption Vr/V a <C 1. Finally, note that 
we make no assumptions about the values of the guide 
field B z and the plasma viscosity. [However, below we 
will see that our assumption of a quasi-stationary recon- 
nection process in a thin current sheet layer results in a 
necessary condition v <C rj (V A /Vn) for the plasma kine- 
matic viscosity v. Refer to Eq. (27) for more details.] 

Now note that several equations that we derived in the 
previous section for the case of the Sweet-Parker recon- 
nection with constant resistivity stay the same in the case 
of anomalous resistivity. Indeed, equation (1) stays valid 
when the resistivity r\ is not constant. Therefore, equa- 
tion (2) also stays valid, except that r] now is the value 
of resistivity at the reconnection layer central point O 
(see Fig. 1), i.e. r\ Q = r/(j z = j Q ,x = 0,y = 0). Equa- 
tions (3) and (8), which result from the Ampere's law 
and the plasma incompressibility respectively, obviously 
stay valid too [39]. At the same time, in the general 
case of anomalous resistivity that we consider in this sec- 
tion, we need to re-derive equations (7) and (9), which 
are the equations of the plasma acceleration and spatial 
homogeneity of the electric field z-componcnt along the 
reconnection layer (i.e. along the y-axis). However, be- 
fore we re-derive these two equations, we would first like 
to derive the equation of magnetic energy conservation. 
Using Eqs. (2), (3) and (8), we immediately obtain 

vB 2 n = e Vo j 2 , cwl, (10) 

where for the purpose of comparison of our theoretical 
results to our numerical simulations we introduce a co- 
efficient e, which is of order of unity. Equation (10) is 
the equation for magnetic energy conservation. The rate 
of the supply of magnetic energy into the reconnection 
layer is equal to the rate of its Ohmic dissipation inside 
the layer. 

Next we derive the equation of plasma acceleration 
along the reconnection layer (i.e. along the y-axis), tak- 
ing into consideration all forces acting on the plasma. 
The MHD equation for the y-component of the plasma 
velocity V y , assuming the quasi-stationarity of the re- 
connection (d/dt = 0) and plasma incompressibility 
(divV = 0), is [24] 

p(V-V)V v = -(d/dy)(P + B 2 x /2 + B 2 y /2) 

+ (B-V)B v + pvV 2 V y , (11) 

where p is the plasma density, v is the plasma kinematic 
viscosity (assumed to be constant) and P is the sum of 
the plasma pressure and the guide field pressure B 2 /2. 
Taking the first order partial derivative d/dy of equa- 
tion (11) at the central point O, we obtain 

pv 2 = -{d 2 P/dy 2 ) +p ]o + pu [W 2 (dV y /dy)] o , (12) 

where we use parameter [3 = (dB x /dy) and the Am- 
pere's law j = (dB y /dx) — (dB x /dy) at point O, and 
we also use the formulas V x = B y = on the y-axis, and 
V y = B x — at point O, which follow from the symme- 
try of the problem with respect to the x- and y-axes. 
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The pressure term (d 2 P/dy 2 ) on the right-hand-side 
of equation (12) can be precisely calculated in analogy 
with the Sweet-Parker derivation of the pressure decline 
along the reconnection layer, which employs the force 
balance condition for the plasma across the reconnec- 
tion layer and leads to equation (4). The viscosity term 
pv [V 2 (dV y j dy)] in equation (12) can be calculated ap- 
proximately by using estimates for the V y velocity deriva- 
tives. In Appendix A we carry out these calculations and 
show that the pressure and viscosity terms are equal to 



(d 2 P/dy 2 ) 



pv [V 2 (dV v /dy)] o 



B m (d 2 B y /dy 2 ) m 

+ o{pv 2 ,0j o ,pvv/S 2 o }, (13) 

-pvv/S 2 , (14) 



where B m (d 2 B y /dy 2 ) m is calculated just outside the re- 
connection layer at point M (see Fig. 1), and expression 
o{pv 2 , 0j o , pvv/8 2 } denotes terms that are small com- 
pared to either pv 2 or 0j o or pvv/8 2 in the limit of a 
slow reconnection rate in a thin reconnection layer, see 
Ref. [40]. Substituting equations (13) and (14) into equa- 
tion (12) and using equation (3) for S D , we obtain 

pv 2 « -B m (d 2 B y /dy 2 ) m +0j o -pvvj 2 /B 2 m , (15) 
where = (dB x /dy) at point O. 

Note that this equation is exact if plasma viscosity can 
be neglected (v — 0). 

Now we use the condition of spatial homogeneity of the 
electric field z-component along the reconnection layer, 
i.e. along the y-axis. We take the second order par- 
tial derivatives d 2 /dy 2 of the left- and right-hand-sides 
of equation (1) at the central point O (note that the 
first order partial derivatives are identically zero). Tak- 
ing into account the symmetry of the problem, so that 
V x = By = on the y-axis, and V y = B x = dj z /dy = 
at point O, we obtain 

- [Vo+jo(dv/dj z ) ] (d 2 Jz /dy 2 ) -j (d 2 r 1 /dy 2 ) 
= 2(dV y /dy) {dB x /dy) = 2v0, 



(16) 



where we use formulas v = (8V y /dy) and = 
(dB x /dy) . Finally, we need to estimate the (d 2 j z /dy 2 ) 
term, which enters the left-hand-side of equation (16). 
This estimation can be done by taking the second or- 
der partial derivative d 2 /dy 2 of equation (3), while keep- 
ing 5 constant because the partial derivative in y is to 
be taken at a constant value x = const — S . In Ap- 
pendix B we give the detailed derivations and find that 
the y-scale of the current j z is about the same as the y- 
scale of the outside magnetic field, i.e. jo 1 (d 2 j z /dy 2 ) w 
B^ l 1 (d 2 B y /dy 2 ) m . However, for the purpose of compari- 
son of our theoretical results to numerical simulations in 
Sec. V, we find it convenient to write 

jo\9 2 j z /dy 2 ) = 1 B m 1 (d 2 B y /dy 2 ) m , 7 « 1, (17) 

where we introduce the coefficient 7, which is of order 
unity. 



Let us take the dimcnsionless coefficients e and 7, 
which enter equations (10) and (17), and are of or- 
der unity, to be exactly unity, e = 1 and 7 = 1. 
Now we have all the equations necessary to determine 
all unknown physical parameters. In particular, using 
Eqs. (10), (15), (16) and (17), we easily obtain the fol- 
lowing approximate algebraic equation for the z-current 
j at the reconnection layer central point O: 



3 + 



jo(dr)/dj. 



Vo 



,, + B m (d 2 v/dy 2 



r, (d 2 B v /dy 2 ) r 



1 + 1 



2 -4 

To3o 



2 Br, 



Vo J V 2 Bi {d 2 B y /dy 2 ) r 



(18) 



where the Alfven velocity Va is defined as Va = B m /y/p 
and r\ = r](j z = j D ,x = 0, y = 0) is the resistivity at 
point O. Given the resistivity function r\ = rj(j z ,x,y), as 
well as the magnetic field B m and its second order deriva- 
tive (d 2 B y /dy 2 ) m outside the reconnection layer, we can 
solve equation (18) for the current j and find the recon- 
nection rate, which is the rate of destruction of magnetic 
flux at point O, equal to — (dA z /dt) . Using Eq. (1), 
we find that the reconnection rate is equal to i] j = E z . 
Note that for the classical reconnection layer that we con- 
sider (see Fig. 1) the right-hand-side of equation (18) is 



positive because 2B m /(d 2 B y /dy 2 ) T 



-L 2 < 0, where 



L is the global scale of the magnetic field outside the 
reconnection layer. Once the current j a is calculated 
by means of equation (18), we can easily calculate all 
other reconnection parameters, using Eqs. (2), (3), (8) 
and (15), 



V R W Wo/Srn « Va, 
V ~ Vo3o/ B m, 

~ jo[(l + ^/Vo)(Vojo/V A B m 
+ (B m /j 2 )(d 2 B y /dy 2 ) r 

So ~ Bra/jo ~ Vr/v. 



< jo, 



(19) 
(20) 

(21) 
(22) 



Equations (18)-(22) are the most general result for mag- 
netic reconnection that we obtain in this paper. Restor- 
ing coefficients e w 1 and 7 w 1, equation (18) becomes 



( 7 + 2e ) + - M d v/ d j*)o + B m (d 2 V /dy 2 



Vo 



1 -4 

%r 



^(d 2 B y /dy 2 ), 
2B m 



Vo J V 2 B? n (d 2 B y /dy 2 ) m 



(23) 



Hereafter we will consider the natural case when 
{drj/dj z ) > and (d 2 r]/dy 2 ) < because plasma con- 
ductivity decreases as the current increases and we are 
interested in fast anomalous reconnection (i.e. faster than 
the Sweet-Parker reconnection). In this case the first, 
second and third terms on the left-hand-side of equa- 
tion (18) are all positive. It is easy to see that the first 
term is related to Sweet-Parker reconnection with con- 
stant resistivity equal to ry , the second term is related 
to fast reconnection associated with the dependence of 
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anomalous resistivity on the current, and the third term 
is related to fast reconnection associated with an ad hoc 
localization of resistivity in space (see the next section for 
details). Also note that if the plasma kinematic viscos- 
ity v is larger than the resistivity rj , then, according to 
equation (18), the current j a and reconnection rate r\ j 
become smaller as v grows, i.e. the reconnection slows 
down for viscous plasmas as one expects. 

We postpone the analysis of equations (18)-(22) until 
Sec. VI. Let us now make an estimate of the half-length 
of the reconnection layer L' (see Fig. 1). Note that L' 
is not needed for the calculation of the reconnection rate 
Vojo by means of equation (18). Nevertheless, we are 
still interested in a rough estimate of L' , in particular, 
because we need to check our assumption that the recon- 
nection layer is thin, i.e. that the condition S Q <C L' is 
satisfied. It is clear that L' can not be much larger than 
the global scale of the magnetic field outside the layer L. 
Therefore we have the condition L' < L. However, L' can 
be much smaller than L, in which case the reconnection 
layer has a pair of the Petschek switch-off MHD shocks 
attached to each end of the layer in the downstream re- 
gions [3, 5, 25], as shown in Fig. 1. In this case L' should 
be calculated as the y-coordinatc of the point on the y- 
axis at which the perpendicular field B x is strong enough 
to support the shocks [5]. Following Kulsrud [5], we use 
the jump condition on the Petschek switch-off shocks to 
obtain [41] 

V R « B x (0, y = L')/y/p « 0L'/y/p = f3L'V A /B m , (24) 

where we use the first-order Taylor expansion for an 
estimate B x (0,y = L') w {dB x /dy) L' = fiL' and, 
as before, Va = B rn j \[p. Now, using Eqs. (19), (21) 
and (24), we can easily find L 1 . Before we write the ex- 
plicit formula for L' , note that the absolute value of the 
(B m / j 2 )(d 2 By I 'dy 2 ) m term in equation (21) is equal or 
smaller than the (1 + 2v /Vo)(Vojo/VAB m ) 2 term. This 
follows from the fact that the left-hand-side of equa- 
tion (18) is equal or greater than unity [see our com- 
ments in the paragraph that follows Eq. (23)]. Therefore, 
the {B m I jt){d 2 By / dy 2 ) m term can be omitted in equa- 
tion (21) for the purpose of estimating V . As a result, we 
obtain the following rough estimates for the reconnection 
layer half-length L' and the velocity of plasma outflow in 
the downstream regions V out : 

l' « (VABl/nofJii + u/no)- 1 

« {Va/v^I + v/^o)- 1 , <5 «L' <L, (25) 

Vout « vL' « V A (1 + vho)- 1 < V A , (26) 

where we use equations (19), (20), (21) and (24). Note 
that the condition V < L is always satisfied because 
2B m /(d 2 By/dy 2 ) m w — L 2 , and the left- and right-hand- 
sides of equation (18) are equal or greater than unity. 
However, the condition that the reconnection layer is 
thin, S <C L', is satisfied only if plasma viscosity is not 
too large, 

v « Vo(V A B m /r ]o j ) w Vo(V A /V R ), (27) 



where, to derive this formula, we use Eqs. (19), (22) 
and (25). In other words, to be able to form a thin 
reconnection layer, the plasma should not be too vis- 
cous. Note that in the case of constant resistivity and 
large viscosity, r\ = const and v 3> Vo, the reconnec- 
tion velocity is V R /V A « (WVa£) 1/2 (W") 1/4 (see [26]) 
and condition (27) reduces to v -C t] (VaL/v ) 2 ^ 3 = 
VaL^o/VaL) 1 / 3 . Finally, if the plasma viscosity is small 
in comparison to the resistivity, v <C Vo , then from equa- 
tion (26) we immediately find an important and well- 
known result that in this case the velocity of the plasma 
outflowing in the downstream regions at the ends of the 
reconnection layer is approximately equal to the Alfven 
velocity, V out w V A (see Fig. 1). 

At the end of this section we would like to discuss sev- 
eral assumptions that we used in our derivations. First, 
the solution of equation (18) is valid only if it gives 
Vojo <C VAB m , which is our assumption of a slow quasi- 
stationary reconnection. Because of equation (19), con- 
dition ri j <C VAB m is equivalent to Vr <C Va, i.e. the 
reconnection velocity, which is the velocity of the incom- 
ing plasma, must be small in comparison to the Alfven 
velocity in the upstream region. Second, the coefficient 
(3 = (dB x /dy) , given by Eq. (21), must be much smaller 
than the current j a = (dB y /dx) — [3 w (dB y /dx) be- 
cause the reconnection layer is assumed to be thin. It is 
easy to see that this condition is satisfied. Indeed, the 
first term in the brackets [...] in equation (21) is much 
smaller than unity because of the upper limit for plasma 
viscosity given by equation (27) . The second term in the 
brackets in equation (21) is also much smaller than unity 
because j is much larger than the electric current out- 
side the reconnection layer due to our assumption of a 
large characteristic Lundquist number of the system. 



IV. SPECIAL CASES OF MAGNETIC 
RECONNECTION 

In this section we focus on three special cases for the re- 
connection rate, which arise when one of the three terms 
on the left-hand-side of equation (18) dominates over 
the other two. We consider the classical Sweet-Parker- 
Petschek reconnection layer shown in Fig. 1 and define 
the global scale of the magnetic field outside the recon- 
nection layer as L = y/—2B m / (d 2 B y /dy 2 ) m (sec [42]). 
In addition, for the purpose of clarity, in this section we 
focus only on resistivity effects and neglect plasma vis- 
cosity, assuming that v <g; r\ . 



A. Sweet-Parker reconnection, r) = const 

First, consider the case when resistivity is constant, 
T](j z ,x,y) = i] = const. In this case only the first 
term on the left-hand-side of equation (18) is nonzero 
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FIG. 2: (Color online) Three top plots: merging of two cylin- 
drical magnetic flux tubes by reconnection (A z is plotted). 
Two bottom-left plots: the configuration of the reconnec- 
tion layer for the Sweet- Parker ( "S-P" ) and Petschek-Kulsrud 
("P-K") cases (j z is plotted). The bottom-right plot: the re- 
connection rate dA z /dt = rj j at the layer central point O, 
normalized by VaoBq- 



and equations (18)-(22) and (25) reduce to 



Z^^VofoL/VABi 
S = V A L/r) > 1, 
v « V A /L, 



1/2 



LSI 



-1/2 



jo ~ (B m /L)S 
V R « V A So 1/2 , 
(3 « (B m /L)So 1/2 , 
V w L, 



(28) 



where we set the plasma kinematic viscosity to zero 
(v = 0), use our definition of the global field scale 
L = y / —2B m /(d 2 B y /dy 2 ) m , introduce the Lundquist 
number S = V A L/r) and assume for our estimates 
that 3 1//4 « 1. The above equations are the Sweet- 
Parker reconnection equations with constant resistivity 
equal to rj Q . Thus, we find that if resistivity is con- 
stant, then the reconnection must be Sweet-Parker and 
not Petschek [5, 22]. We discuss this important result in 
Sec. VI. Note that in this section, contrary to Sec. II, 
we do not assume the Sweet-Parker geometrical configu- 
ration for the reconnection layer, but derive it together 
with the reconnection rate from our general equations of 
the previous section. A typical configuration of the re- 
connection layer in the case of Sweet-Parker reconnection 
is shown on the left-bottom plot in Fig. 2. This plot is 
marked by letters "S-P" . 



B. Petschek-Kulsrud reconnection, r\ = ri(j z ) 

Now let us consider the case when resistivity is anoma- 
lous and is a monotonically increasing function of the 
electric current only, r\ = rj(j z ). Let us further assume 
that this dependence of resistivity on the current is very 
strong, so that {jo/Vo)(dr]/dj z ) 1. In this case the sec- 
ond term on the left-hand-side of equation (18) is domi- 



nant, and equations (18) and (25) reduce to 

,3-3r2/T/2 D 4 _ 



(dv/dj z )o 
Vr 
V a 

L' 



Vojo _ <5c 

V A B m ~ L 
VaBI 



B„ 



drj 



Vojo 



L 



V A L 2 \d Jz J o 
(V A j Lyr, B m y/* 



1/3 



L 



jo ( dr) 
Vo \dj z 



-1/2 



< L. 



(29) 



(30) 



Here again we set v — 0, and we use the formula 
L 2 = -2B m /(d 2 B y /dy 2 ) m and Eqs. (19) and (22). From 
equation (30) we see that the half-length of the recon- 
nection layer L' is much less than the global field scale 
L in the case of a strong dependence of resistivity on 
the current ((j /Vo){dr]/dj z ) 3> 1). This means that 
in this case the geometry of the reconnection layer is 
Petschek with a pair of shocks attached to each end of 
the layer in the downstream regions, as shown in Fig. 1. 
Equation (29) for the reconnection rate was first ana- 
lytically derived by Kulsrud [5] for a special case when 
dri/dj z = const = 77*/ j c . In the Petschek-Kulsrud recon- 
nection case the reconnection rate, given by Eq. (29), can 
be considerably faster than the Sweet-Parker reconnec- 
tion rate. This fast reconnection has been previously ob- 
served in many numerical simulations done with anoma- 
lous resistivity 77 = T)(j z ) (e.g. see [9, 10, 27]). The typ- 
ical configuration of the reconnection layer in the case 
of Petschek-Kulsrud reconnection is shown in the left- 
bottom plot in Fig. 2, this plot is marked by letters "P- 
K". 



C. Spatially localized reconnection, 77 = ri(x,y) 

Finally let us consider the special case when resistiv- 
ity is given by 77 = r)(x,y) and it is spatially localized 
around the central point O of the reconnection layer 



-2i lo /{d 2 r 1 /dy 2 ) 



= I 2 

— s 



(see Fig. 1), so that 

—2B m /(d 2 B y /dy 2 ) m . In other words, we assume that 
resistivity is anomalous and is localized on scale l v that 
is much smaller than the global field scale L. In this case 
the third term on the left-hand-side of equation (18) is 
dominant, and equations (18)-(22) and (25) reduce to 



L/l n « 7 h j 2 L/VAB 2 m 
Si = V A l v /r) > 1, 

V W V A /l v , 
So 



jo ~ (Bjn/l^S^ ^ , 



-1/2 



-1/2 



Vr W V A S i 
13 « (B m /l v )S l 
V w l v < L, 



-1/2 



(31) 



where we again set v — 0. The above equations are the 
same as the Sweet-Parker equations (28) with the global 
field scale L replaced by the resistivity scale l v <C L. Note 
that, when resistivity is localized, the reconnection rate 
becomes faster than the Sweet-Parker rate by a factor 
yjL/l n ^> 1 and the geometry of the reconnection layer 
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is Petschek with a pair of shocks attached to each end 
of the layer (see Fig. 1). These results are in agreement 
with many previous numerical simulations of reconnec- 
tion with spatially localized resistivity [7, 8, 11, 28]. 

V. NUMERICAL SIMULATIONS OF 
MAGNETIC RECONNECTION 

In this section we present the results of our numerical 
simulations of unforced reconnection of two cylindrical 
magnetic flux tubes. These simulations are not intended 
as a check or a proof of our theoretical results for mag- 
netic reconnection. Our equations (18)-(22) have been 
derived analytically and are very general. A compre- 
hensive testing of them would require extensive compu- 
tational work, which is beyond the scope of this paper. 
Instead, we present our simulations as a demonstration 
of our reconnection model predictions. 

Following Kulsrud [5] , we assume that plasma resistiv- 
ity is given by the following piecewise linear function of 
the z-component of the electric current: 

V(jz) =Vs+V* max{0, j z - j c }/jc, (32) 

where r} s is the Spitzer resistivity, which is assumed to 
be very small, 77* is the anomalous resistivity parame- 
ter and j c is the critical current parameter. The Kul- 
srud model's prediction for the reconnection rate in the 
case (jo/Vo){dv/dj z ) = joV*/Vojc > 1, which is given by 
equation (29), has already been checked and confirmed 
numerically by Breslau and Jardin [27]. Here we simulate 
reconnection with anomalous resistivity given by equa- 
tion (32) for a different computational setup, a higher 
Lundquist number and a wider range of parameters n+ 
and j c (without the restriction joV+fVojc 3> 1)- Our in- 
tent is to see how our general formula (18) for the recon- 
nection rate works in this case. 

We consider an unforced reconnection of two cylindri- 
cal magnetic flux tubes with the initial z-component of 
the magnetic field vector potential equal to 

A z (x, y) = A [cxp (-r 2 + /2R 2 ) + cxp (-r 2 _/2R 2 )] , 
A = B RoV^, r± = (xTd) 2 +y 2 , (33) 

see the top-left plot in Fig. 2. This convenient compu- 
tational setup was suggested to us by Mikic and Vain- 
shtein [29]. We choose the parameters in equation (33) as 
Ro = 1 (the global scale of the field is unity), d = 2.7162 
(5% of initially reconnected flux) and B = 1 (the max- 
imal initial field is unity). Thus, A n = \fe. In addition, 
for convenience we choose the plasma density p = 1, 
so that the typical Alfven velocity and time are unity, 
Vao = -Bo/VPo = 1 and t ao = Rq/Vao = 1- The guide 
field is chosen to be zero, B z = 0, and the initial plasma 
velocities are zero. The initial gas pressure P is cho- 
sen in such way that each of the two cylindrical magnetic 
flux tubes (33) would initially be in complete equilibrium, 
P+B 2 /2+B 2 /2 = const, if there were no magnetic forces 



from the other tube. The plasma kinematic viscosity is 
chosen to be equal to the Spitzer resistivity, v = r\ s . The 
boundary conditions are placed at x, y = ±25, which are 
virtually at infinity (the magnetic vector potential (33) 
drops to less than 10~ 100 at the boundaries). Because 
of the symmetry of the problem, in the case of a quasi- 
stationary reconnection considered here, it is enough to 
run simulations only in the upper-right quadrangle of the 
full computational box. We use the FLASH code for 
our simulations. This is a compressible adaptive-mesh- 
refincment (AMR) code written and supported at the 
ASC Center of the University of Chicago. (For a compre- 
hensive description of the FLASH code see [30, 31]). The 
MHD module of the code uses central finite differences to 
properly resolve all resistive and viscous scales. Compar- 
ing numerical results obtained by simulations done with 
a compressible code to our approximate theoretical for- 
mulas derived for incompressible fluids is fine in a case 
of a very high Lundquist number. This is because in 
this case the incompressibility condition is a very good 
approximation in a reconnection layer even for compress- 
ible plasmas [23]. Indeed, in our simulations the plasma 
density varies by no more than 15% inside the reconnec- 
tion layer. The biggest advantage of the FLASH code for 
our purposes is that it is an already existing, well tested 
code with the AMR feature, which allows us to place the 
boundary conditions at infinity. The size of the smallest 
elementary grid cell in our two-dimensional simulations 
was chosen typically to vary from 25/2 15 = 0.0007629 
to 25/2 13 = 0.003052, which was sufficient to resolve the 
resistive reconnection layer. 

The two cylindrical magnetic flux tubes, initially set 
up according to equation (33), attract each other (in a 
similar way as two wires with colliniar currents do). As 
a result, as time goes on, the tubes move toward each 
other, form a thin reconnection layer along the y-axis 
and eventually completely merge together by reconnec- 
tion. This merging process is displayed in the three top 
plots in Fig. 2, which show the field vector potential A z 
in a central region of the full computational box. The 
two bottom-left plots in Fig. 2 show the electric cur- 
rent j z in a central region that includes the reconnection 
layer. These plots clearly show the reconnection layer 
configuration which is formed during the reconnection 
process in the cases of the Sweet-Parker and Petschek- 
Kulsrud reconnection (refer to Sees. IV A and IV B). The 
bottom-right plot in Fig. 2 demonstrates the functional 
dependence on time typical of the normalized reconnec- 
tion rate (1 /VaoBo)cIA z / dt = r] j /V A oBo at the recon- 
nection layer central point O. (This point is shown in 
Fig. 1). Next we compare the maximal (peak) reconnec- 
tion rate observed in the numerical simulations with the 
theoretical rate predicted by our reconnection model. 

When resistivity is given by equation (32) and v = 
r] s , our theoretical formula (18) for the reconnection rate 
reduces to 

3 I Vlk ~n I „ /„ g o 6 12 P (34) 

in ~ { L + r, s/Vo) y 2 R4 R6 R2 , {^) 
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FIG. 3: Three left plots: the reconnection rate (left-top), co- 
efficient e (left-central) and coefficient 7 (left-bottom) as func- 
tions of 77* at fixed 77^ = 0.0002 and j c — f7.f00. The crosses 
and boxes show the results of our simulations and theory re- 
spectively. The solid/dashed horizontal lines correspond to 
77* = and simulations/theory. The inclined dotted line given 

1/3 

in the reconnection rate plot shows the oc rjj scaling. Point 
M (see Fig. 1) is chosen to satisfy (r)j z ) m = (l/3)r} jo- Three 
right plots: the same as the three left plots except point M is 
chosen to satisfy (ijj z ) m = (l/W)riojo- 



where, as explained above, in our numerical simula- 
tions we choose B = 1, R = 1, Vao = 1 and 
po = 1, and by definit ion the global field scale is L = 
y/—2B m /(d 2 B y /dy 2 ) m . We also assume that j a > j c , 
which is the case that we consider in our numerical simu- 
lations. As we can see from equation (34), the theoretical 
reconnection rate r\ j is rather sensitive to B m and L, 
which are the strength and scale of the magnetic field 
at point M outside the reconnection layer (see Fig. 1). 
Therefore, in order to compare the theoretical results 
with the results of our simulations, we need to accurately 
calculate the B m and L observed in the simulations. As 
a result, the choice of the exact position of point M, at 
which B m and L are calculated, is important. First, in 
our simulations we choose the point M to be the point 
on the x-axis at which the observed resistivity term r\j z 
is three times smaller than that observed at the central 
point O, i.e. (r]j z ) m = (l/3)r) j . The three plots on the 
left in Fig. 3 demonstrate our results for this choice. The 
top plot on the left shows the reconnection rate. The 
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FIG. 4: These plots are the same as those in Fig. 3, except 
here j c = 4.135. 



crosses are a log-log plot of the maximal (peak) recon- 
nection rate observed in our simulations as a function 
of parameter 77* for fixed r] s = 0.0002 and j c = 17.100. 
The boxes are the theoretical reconnection rate, which 
is given by equation (34) with the appropriate values of 
B m , L and density p observed in the simulations. The 
solid horizontal line (simulations) and the dashed hori- 
zontal line (theory) correspond to the 77* = case. The 

1/3 

inclined dotted line demonstrates the oc r\J scaling [re- 
fer to Eq. (29)]. The crosses and the boxes do not follow 
the oc r]l^ 3 scaling for large values of 77* simply because 
in this case the reconnection rate becomes relatively fast 
and the magnetic field B m outside of the reconnection 
layer is not piled up as much as in the case when 77* is 
small and the reconnection rate is relatively slow. As 
a result, our rate curves flatten at large values of 77*. 
The central and the bottom plots on the left in Fig. 3 
show the observed values of coefficients e and 7, which 
are directly calculated from the simulation data by using 
equations (10) and (17). As we can see, e and 7 are of 
order of unity, as one expects. The three plots on the 
right in Fig. 3 are the same as the three plots on the 
left except point M is now chosen as the point on the 
x-axis at which {f]j z )m = (l/10)r]ojo- Comparing the 
plots on the left and the plots on the right, we see that 
the choice of point M is indeed important. Note that for 
the choice (r]j z ) m = (l/3)r] j for the position of point 
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M the half-thickness of the reconnection layer S , defined 
as the abscissa of point M, increases from 0.011 to 0.050 
as increases from zero to its maximal value shown 
on the plots in Fig. 3. At the same time, for the choice 
(V3z)m = (l/10)?7ojo for the position of point M S a ranges 
from 0.018 to 0.056, which is noticeably larger. Perhaps 
simulations of forced reconnection with a strict control 
of position of point M together with control of the out- 
side field B m and its scale L, could be better suited for 
comparison to our theoretical model. Such simulations 
are beyond the scope of this paper. However, see more 
discussion of forced reconnection in the next section. Fi- 
nally, Figure 4 shows the same results as Fig. 3, except 
the former has plots for a smaller value of the critical 
current, j c — 4.135. 

We believe that the results presented in Figs. 3 and 4 
generally confirm our theoretical model. In particular, in 
all cases the theoretical reconnection rates and the rates 
observed in the simulations do not differ by more than 
33%. The observed relatively small discrepancy between 
the theoretical and simulated rates is mainly due to the 
coefficient e not being precisely constant in the simula- 
tions, while the variations of coefficient 7 are somewhat 
less important [see Eq. (23) for the theoretical reconnec- 
tion rate]. This discrepancy can be due to two causes. 
First, our theoretical model is general, but approximate, 
and second, the plasma is compressible in the simula- 
tions, while it is assumed to be incompressible in the 
theoretical model. 



VI. DISCUSSION AND CONCLUSIONS 

Let us summarize our main results. In this paper we 
take a new theoretical approach to the calculation of the 
rate of quasi-stationary magnetic reconnection. Our ap- 
proach is based on analytical derivations of the recon- 
nection rate from the resistive MHD equations in a small 
region of space that is localized about the center of a 
thin reconnection layer and has its size equal to the layer 
thickness. Our local-equations approach turns out to be 
feasible and insightful. It allows us to consider mag- 
netic reconnection with an arbitrary anomalous resistiv- 
ity and to calculate the reconnection rate for this general 
case [see Eq. (18)]. We find the interesting and impor- 
tant result that if plasma is incompressible and reconnec- 
tion is quasi-stationary, then the reconnection rate is de- 
termined by the anomalous resistivity function rj(j, x, y) 
and by the strength and structure of the magnetic field 
just outside of the reconnection layer (i.e. at point M in 
Fig. 1). Thus, we find that the global magnetic field and 
its configuration are not directly relevant for the purpose 
of calculation of a quasi-stationary reconnection rate, al- 
though, of course, the local magnetic field outside the 
reconnection layer depends on the global field. 

One of the major results of this paper is that in the 
case of constant resistivity, i] = const = rj , the mag- 
netic reconnection rate is the slow Sweet-Parker recon- 



nection rate and not the fast Petschek reconnection rate 
[refer to Eqs. (28)]. This result agrees with numerical 
simulations and at the same time contradicts the re- 
sult of the original Petschek theoretical model. In the 
framework of our theoretical approach, based on local 
equations, the reason for this contradiction can be un- 
derstood as follows: In the Petschek model our param- 
eter v = (dV y /dy) = —(dV x /dx) , which is equal to 
the first order partial derivatives of the incompressible 
plasma velocities at the reconnection layer center, is ba- 
sically treated as a free parameter. This is because in 
the Petschek model v can be estimated as the ratio of 
the plasma outflow velocity (equal to the Alfven veloc- 
ity for viscosity-free plasma) and the reconnection layer 
length, v = (dV y /dy) w Va/L', and the layer length 
L' is treated free parameter by Petschek. In his 
model L' is taken to be equal to the minimal possible 
value, that the Petschek shocks do not seriously perturb 
the magnetic field in the upstream region. This is L' w 
(L/So) (In S ) 2 , where S = VaL/t] is the Lundquist 
number [3, 5]. In this case v w {V^/r] ) (In S Q )~ 2 and, ac- 
cording to equations (19) and (20), the reconnection ve- 
locity in this case is equal to Vr w V^lnS^) -1 , which is 
the Petschek result. On the other hand, in our theoretical 
model the parameter v is not treated as a free parameter. 
In fact, our three physical parameters v = {dV y /dy) 0l 
[3 = (dB x /dy) and j D = j z (x = 0,y = 0) are con- 
nected to each other and must be calculated from equa- 
tions (10), (15) and (16). Let us discuss the meaning of 
these equations. Equation (10) is the equation of mag- 
netic energy conservation. It says that the rate of sup- 
ply of the magnetic field energy into the reconnection 
layer vB^ must be equal to the rate of the resistive dis- 
sipation of this energy r) j 2 . In the Petschek model v 
and, accordingly, the rate of the magnetic energy supply 
vB^ are basically prescribed by hand (resulting in an ad 
hoc fast reconnection), while in our model they are self- 
consistently calculated from the MHD equations. Equa- 
tion (15) is the equation of plasma acceleration along 
the reconnection layer. It says that the magnetic tension 
force (3j must be large enough in order to be able to push 
out all the plasma along the layer that is supplied into 
the layer. Finally, equation (16) is the equation of spa- 
tial homogeneity of the electric field z-component along 
the reconnection layer. This equation sets an upper limit 
on the product vj3 in the case of a quasi-stationary re- 
connection and is directly related to the calculations and 
arguments given by Kulsrud [5] in the framework of the 
global-equations theoretical approach, sec Rcf. [43]. As a 
result, none of parameters v, [3 and j a can be treated as 
free parameters, and all of them must be self-consistently 
determined from the MHD equations. 

It is very instructive to briefly examine our results from 
the two distinct points of view in connection with the 
reconnection problem, which are expressed in numerous 
papers on computer simulation of magnetic reconnection. 
These two points of view are: unforced (free) magnetic 
reconnection and forced magnetic reconnection. In the 
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case of unforced reconnection, one should solve our main 
equation (18) for the current j a at the reconnection layer 
center and then solve for the reconnection velocity Vr 
by using equation (19). The solution for j D and Vr will 
depend on B m , which is the strength of the magnetic 
field outside the reconnection layer and which enters the 
right-hand-side of equation (18). On the other hand, in 
the case of forced magnetic reconnection the reconnec- 
tion velocity Vr is prescribed and fixed. In this case the 
magnetic field outside the reconnection layer B m should 
be treated as an unknown quantity, and equations (18) 
and (19) should be solved together in order to find the 
correct quasi-stationary values of j and B m . In other 
words, in the forced reconnection case an initially weak 
outside magnetic field B m gets piled up to higher and 
higher values until the resulting current j a in the recon- 
nection layer becomes large enough to be able to exactly 
match the prescribed reconnection velocity Vr and to be 
able to reconnect all magnetic flux and magnetic energy, 
which are supplied into the reconnection region in the 
quasi-stationary reconnection regime. 

Finally, a couple of words about plasma viscosity and 
guide field and their effect on magnetic reconnection. 
First, according to our equations (18) and (19), in the 
case when the resistivity is constant, rj = const = 
7] , and the plasma viscosity is much larger than re- 
sistivity, v » 7] , the reconnection velocity becomes 
Vr/Va ~ S 1 ^ 2 {rj /v) 1 / i , which is \fvjr\ times smaller 
than the Sweet-Parker reconnection velocity given by for- 
mula (28), see [26]. Thus, we see that the reconnection 
rate becomes smaller when the plasma viscosity becomes 
large. However, note that in many astrophysical and 
laboratory applications plasmas are very hot and highly 
rarefied. Under these conditions the ion gyro-radius be- 
comes much shorter than the ion mean-free-path, and the 
plasma becomes strongly magnetized. As a result, the 
plasma viscosity becomes the Braginskii viscosity, which 
is dominated by magnetized ions [32]. In this case in 
all our equations above the isotropic viscosity v, which 
is proportional to the ion mean-free-path, should be re- 
placed by the Braginskii perpendicular viscosity, which is 
proportional to the ion gyro-radius and is much smaller 
than the perpendicular viscosity in a strongly magne- 
tized plasma. Second, according to our results, in two- 
and-a-half dimensional geometry the guide field B z has 
no effect on the quasi-stationary reconnection rate. In- 
deed, in our derivations the guide field appears only as 
magnetic pressure B z /2 term in addition to the plasma 
pressure. The combined pressure P enters equations (11) 
and (12) of plasma acceleration along the reconnection 
layer and the value of spatial derivative of P is given 
by equation (13), which docs not involve B z . Thus, the 
guide field gets eliminated and does not enter into our 
final equation (18) for the reconnection rate. However, 
if one assumes that the anomalous resistivity rj depends 
on x- and y-componcnts of the current j x = dB z /dy and 
j y = —dB z /dx in addition to its dependence on the z- 
componcnt of the current j z , then the reconnection rate 



will depend on the guide field B z . 

In this paper we consider quasi-stationary magnetic 
reconnection in a thin reconnection layer. We leave a 
study of tearing modes instability in a reconnection layer 
and non-quasi-stationary reconnection for a future paper. 
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APPENDIX A: DERIVATION OF 
EQUATIONS (13) AND (14) 

Below, for brevity, we assume that spatial derivatives 
are to be taken with respect to all indexes that are listed 
after the comma signs in the subscripts, e.g. V x yy = 
d 2 V x /dy 2 . 

We derive equation (13) first. The derivation is some- 
what analogous to the Sweet-Parker derivation of the 
pressure decrease along the reconnection layer, which 
leads to equation (4). Namely, to find the pressure de- 
crease and the pressure second derivative along the layer, 
we integrate the pressure gradient vector along the con- 
tour O— »M— >M— »0 shown in Fig. 1 and use the force 
balance condition for the plasma across the reconnection 
layer. Now we carry out these calculations in a mathe- 
matically precise way. 

For infinitcsimally small values of the y-coordinate, 
taking into account the symmetry of the reconnection 
layer with respect to the x- and y-axes and plasma in- 
compressibility, we use the following Taylor expansions 
in y for the plasma velocity V(ir, y) and for the magnetic 
field B(x,y): 

V x = Ui 0) (x) + tf/2)V^ y {x) + (y 4 /24)ui°U(x), 
B x = yB x °} y {x) + ( y y6)B x % y (x), 
V y = -yV^(x) + ( y y6)V y { %(x), 
B y = B y °\x) + ( y y2)B y a) yy (x), 
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where the variables with the superscripts ^ are taken at 
y = and depend only on x. Assuming quasi-stationarity 
of reconnection (d/dt = 0) and plasma incompressibility, 
the MHD equation for the plasma velocity V is 

VP + V(£? 2 + By)/2 = -p(VV)V + (BV)B 

+/^V 2 V, (A2) 

where P is the sum of the plasma pressure and the z- 
componcnt field magnetic pressure B 2 /2. Let us cal- 
culate line integrals of the left- and right-hand-sides of 
equation (A2) along the contour O— >M— >M— >0 shown 
in Fig. 1. First, the line integral of the left-hand-side is 
obviously 

P(0,y)-P(0,0) + (1/2)^(0, y) 

= (y 2 /2)[P iVV (0,0)+n (A3) 

where y is the y-coordinate of points M and O, and we 
use the formulas B y = on the y-axis, B x = at point 
O and B x (0 7 y) = /3y for small y [see the definition of /3 
in Eq. (15)]. Second, using expansion formulas (Al), we 
calculate the following line integrals along the contour 
O— »M— >M— »0, up to the second order in y: 

|[ P (vv)v]di = ^([vW] 2 -^[tl + t]] 

r M 

+ 2 J VMvg^dx), (A4) 
|[(BV)B]dl = ^(BWBH+/3 2 +BWjW 

/ I V x.yyy 
JO 

+ B i>yV d *)> ( A5 ) 

J [^v 2 v]di = (2v^] y v^l - yw w (0) 

-J V *%yy dx )> ( A6 ) 

where the variables with the superscripts ( m ) are calcu- 
lated at point M (see Fig. 1). 

Next we estimate the terms on the right-hand-sides 
of Eqs. (A4)-(A6). Recall the notations V and 6 a for 
the half-length and half-thickness of the reconnection 
layer (see Fig. 1), and that S /L' <C 1 (which is our 
assumption of a thin reconnection layer). The z-current 
inside the reconnection layer is approximately equal to 
jo ~ B m /5 Q [see Eq. (3)], where B m = B y m ^ is the 
field at point M, while the z-current outside the layer 
is ji m ^ < B m /L'. The last two terms on the right-hand- 
side of equation (A5) can be estimated as B x ™^ j^ ~ 
0B m /L' ~ 0jo(S o /L') « (3 j Q and $ [B^ + 

Bx 0) y B y 0) yy }dx ~ 5 B m (i/L 12 ~ (3j (6 /L') 2 « 0j o . Next, 
from equation (1) with the resistivity term dropped we 



see that outside the reconnection layer the typical scale 
of the plasma velocity V is about the same as the typ- 
ical scale of the magnetic field B and, therefore, can 
not be smaller than L' . Thus, the estimates of the 
three terms on the right-hand-side of equation (A4) are 

[V^f ~ 14 (m) [Fi™2 + V™] ~ V 2 /L' 2 ~ v 2 (8 /L') 2 « 
v 2 [see Eq. (8)] and J™ V^V^y dx ~ 5 vV R /L' 2 ~ 
v 2 {8 /L') 2 <C v 2 . The estimates of the four terms on the 
right-hand-side of equation (A6) are V y ™ yy ~ V x ™l x ~ 
1V£' 3 ~ ^o/i' 3 « V y ( % y (0) ~ t,/£' 2 « W /<5 2 

and vj ^ rfx ~ SoVr/L" 1 ~ W <5 2 /L' 4 « W /<5 2 . Note 
that here we use Z/ for estimation of y-derivatives. In 
fact, using the global scale L of the magnetic field outside 
the reconnection layer would have been more appropriate 
for some of the estimations (as shown in Appendix B). 
However, using 1/L' > 1/L for the upper estimates of 
the d/dy derivatives is perfectly fine for the purposes in 
this appendix. 

Next, calculating the line integral of the right-hand- 
side of Eq. (A2) by using formulas (A4)-(A6) and our esti- 
mates made in the previous paragraph, and then equating 
the result to the right-hand-side of Eq. (A3), we easily ob- 
tain equation (13). Note that the first term on the right- 
hand side of equation (13) can also be written in terms of 
the second derivative of the magnetic pressure at point 
M: B m (d 2 B y /dy 2 ) m = (d 2 /dy 2 )(B 2 /2 + B 2 /2) m - (3 2 = 
{8 2 1 dy 2 ){B 2 y /2 + B 2 x /2) m + o{(i ]o } (note that (3 < j a be- 
cause the reconnection layer is thin). Therefore, as noted 
by Zweibel [33], equation (15) is similar to Bernoulli's 
equation. 

Now we derive equation (14), which gives an approxi- 
mate estimate of the viscosity term pv [V 2 (dV y / dy)] o in 
equation (12). We make this estimate as follows: Note 
that V yt y X = on the y-axis because of the symmetry 
of the problem relative to this axis. Therefore, from the 
second order Taylor expansion of V ytV (x, 0) in x, we ob- 
tain an approximate formula V y _ yxx {0, 0) ~ \Vy,y{$oi 0) — 
Vy, y {°>°)}/ S o ~ -Vy, y {0,0)/5l = -v/5 2 , where we take 
into account that V y , v (S o ,0) <~ Vr/L' w vS /L' <C v. 
We can rewrite this approximate formula for V y , yxx (0,0) 
as the following exact formula: V ytVXX (0, 0) = —Cv/5 2 , 
where C is an unknown coefficient of order unity. Our 
numerical simulations of reconnection with constant re- 
sistivity show that C is indeed about unity if S is es- 
timated by equation (3). Thus we immediately find 
that pv [V 2 (dV v /dy)] o = pv[V y , yxx (0,0) + V y . yyy (0,0)} « 
pvV Vtyxx (0, 0) w —pvv/5 2 , which is equation (14). Here 
we also use V y:Vyy (0,0) <~ v/L' 2 <c v/5 2 . 



APPENDIX B: DERIVATION OF EQUATION (17) 

Here as in Appendix A, we assume that spatial deriva- 
tives are to be taken with respect to all indexes that 
are listed after the comma signs in the subscripts, 
e.g. B XjXy = d 2 B x /dxdy. We derive equation (17) in 
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two steps. 

First, we estimate B x yyy at points O and M (see 
Fig. 1), since we will need these estimates below. Con- 
sider the formula B x , x — —B y ^ y , which represents the 
fact that the magnetic field is divergence-free. Take the 
d/dy and d 3 /dy 3 derivatives of this formula and integrate 
the resulting equations along the interval OM shown in 
Fig. 1. We obtain 



The integral Jq B x ^ yyy dx can be estimated by using for- 
mula (B3). The integral of j z . yy can be estimated as 



So 3z, yy dx ~ S jz°yy As a result, we obtain 



u oJ z ,yy ~ D y,yy ^ °° ^2 ° r o D « 



' L 2 ' 



y,yy 



r(™) 
y,yy 



(B5) 



r("») _ r(°) 
x,yyy x,yyy 



I B x yx dx / L 
Jo Jo 



By,yy dx 



-S /? (m) 
u °"y-,yy 



la' 



-E$x,yyyx dx 



/ B y 
Jo 



a mm 



° y,yyyy> 



where the variables with the superscripts (°) and ( m ) are 
taken at points O and M respectively and (3 — B x °]j [see 
Eq. (15)]. Note that S is the half-thickness of the re- 
connection layer, equal to the abscissa of point M (see 
Fig. 1). In making the estimates of the integrals on 
the right-hand-sides of Eqs. (Bl) and (B2), we take into 



(o) 



B 



(o) 



because B y = 



account that B Vt 

on the y-axis. Now, using Eq. (Bl), we estimate that 
B x "^ w P — 6 B y ™\. Let L be the global scale of the 
magnetic field outside the reconnection layer. Then, since 
point M is located outside the reconnection layer, we have 
B^yyyy ~ ±B y ™i/L 2 and B x m y ] yy ~ ±B x m y / L 2 . Next, 
using these estimates, the above estimate for B x m y and 
Eq. (B2), we obtain the following formula: 

Bii m - B x ^ yy ~ ±/3/L 2 ± (5 /L 2 )B^> v . (B3) 

Second, we estimate j z ,yy at point O. Consider Am- 
pere's law formula j z = B VyX — B x , y . We take the d 2 /dy 2 
derivative of this equation and integrate the result along 
the interval OM shown in Fig. 1. We find that 

/■M i-M i-M 

I Jz.yy dx / By^yy X dx / B x yyy dx 
Jo Jo Jo 



= n( m ) _ 
y-yy 



M 



B Xj yyy dx . 



where we use S Q <C L and 5 f3/L 2 = (5 j /L 2 )(f3/j ) 



(Bl) (B m /L 2 )(P/ Jo ) w ±B v m v l{(3/j ) « \BZ'y\, see Eq. (3) 
and note that (3 j because reconnection layer is 
dx thin. Finally, substituting an estimate S w B m /j a into 
formula (B5), we obtain equation (17) with coefficient 

(B2) 7-1- 

As suggested by Uzdensky [34], there exists a nice 
graphical interpretation of the fact that the y-scale of the 
current j z is about the same as the scale of the outside 
magnetic field, i.e. j /(d 2 j z /dy 2 ) « B m /(d 2 B y /dy 2 ) m 
and that equation (17) holds. There can be two different 
cases of the reconnection layer geometry. First, the half- 
length of the reconnection layer V can be approximately 
equal to the global scale L of the outside field. In this 
case the reconnection is Sweet-Parker and V w L is the 
only available scale in the y-direction. Therefore, in this 
case j /\d 2 j z /dy 2 \ w B m /\d 2 B y /dy 2 \ m w L 2 w L' 2 . In 
the second case the reconnection layer half-length is much 
smaller than the global scale, L' L, and the reconnec- 
tion is fast (relative to the Sweet-Parker reconnection). 
In this case, the z-current j z (0,y) on the y-axis drops 
abruptly, as the y-coordinate passes value V and point 
(0, y) moves from the region inside the reconnection layer 
to the region of the outflowing plasma that is located be- 
tween the Petschek shocks (see Fig. 1). However, the z- 
current j z stays large inside the shocks. In other words, 
j z is a smooth function (on the global scale L) along the 
lines that lie inside the reconnection layer and extend 
into the shock separatrices. Thus, in this case, despite 
L'<L, the y-scale of the z-current at the reconnection 
layer central point O is still L. This graphical interpreta- 
tion is well demonstrated by the bottom-left plot of the 
current for the Petschck-Kulsrud ("P-K") reconnection 
(B4) case in Fig. 2. 
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